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The spectrum of vacuum fluctuations in the Milne space (i.e. the t — p coordi¬ 
nate system) is an important ingredient in the thermalization studies in relativistic 
heavy ion collisions. In this paper, the Schrodinger functional for the gauge theory 
perturbative vacuum is derived for the Milne space. The Wigner-transform of the 
corresponding vacuum density functional is also found together with the propagators. 
We finally identify the fluctuation spectrum in vacuum, and show the equivalence 
between the present approach and the symplectic product based method [1 


Q. 


INTRODUCTION 


When two heavy ions collide at ultra-relativistic energies, the collision produces a system 


that will even 


initial stage 


ua 


ly evolve into quark-gluon plasma. In the CGC-Glasma picture of the 


13| . each nucleus is initially composed of two parts; a thin disk made up of 


large x partons acting as a static colour source and the strong static gluon field generated 
by the disk. Gluon quanta in this gluon field are the small x partons in the nuclear parton 
distribution. When the nuclei collide, the large x partons mostly pass through each other 
almost unscathed (except occasionally producing hard jets) but the gluon helds from each 
nuclei start to interact strongly with each other. 

In the lab frame, the two Lorentz contracted nuclei start to overlap at t ~ z ~ 0. At any 
t > 0, the system is then composed of the two disks located at ^ ~ and the gluon field 
stretched between them. In a very simplified picture, one could think of the gluonic system 
as a uniform ‘string’ that is being stretched between two ends pulling away with the speed 
of light. In this ‘string’ picture, a segment of string at position 2 ; at time t will have the 
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speed V = z/t. Since the length of the string is linearly proportional to time t, this speed 
will remain constant for any given segment of the string. This segment therefore can be 
uniquely labelled by the space-time rapidity 

rj = tanh~^ {z /1) ( 1 ) 

and its local time is the proper time 

T = t/'y = — z‘^ ( 2 ) 


which are nothing but the Milne coordinates. Therefore, this coordinate system is the most 
natural one to describe the evolution of such systems. Note that at r = 0 the spacetime 
rapidity rj is not well dehned. Accordingly, theory of the created matter must be formulated 
in this space with the restriction that r is strictly non-zero and positive dehnite. 

The usual perturbation theory in Minkowski space is concerned about calculating the 
scattering amplitude M. = (out|T|in) where T is the scattering operator, the |in) state is 
dehned at the remote past t = —oo, and the |out) state is dehned at the remote future 
t = -l-cxD. Hence, the problem being solved is not an initial value problem but a boundary 
value problem. In contrast, the quantum held theory in the Milne space is most naturally 
formulated as an initial value problem since r is restricted to be positive, and also r = O’*" 
(the trajectory of the colour sources) represents an actual physical boundary in any ultra- 
relativistic heavy ion collisions. 

The initial value problem in quantum mechanics is most naturally formulated in terms 
of the expectation value (0(r)) = (in|t/l(r)Of/(r)|in) where U{t) is the time evolution 
operator. Computation of the expectation value requires the path integral on the closed 
time path (CPT)[m, 15 |d Explicitly, 






(6(r))= [dct>f] / 1201 / (3) 


where 


Po[0?,02] = (0i| in) (in| 02> (4) 

is the matrix element of the initial density operator po = I in) (in |. The details of O[0i,02] 
depends on the time-ordering structure of the operators in O. The held 0i lives on the 

^ This formalism is variously known as the in-in formalism, Schwinger-Keldysh formalism, and also Keldysh- 
Schwinger formalism. 
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forward going time line and the field 02 lives on the backward going time line. Here tq > 0 
is the initial time, and all the quantities with a 0 subscript are evaluated at Tq. The two 
helds 01,2 share the same boundary value 0/ at the hnal time which is then traced. Here, 
the functional integral measure P0 represents integrating over a function of both space and 
time and [d0-^] represents integrating over a function of space at a hxed time (in this case 
t = tf). 


Throughout this paper, we denote the position and the position integral in the (longitu¬ 
dinal) Milne space with a tilde. That is, x = (x_l, rj) and d^x = d^xxdr]. The subscript T 
denotes only the x and y components of a vector quantity. 

To compute (0(r)) in Eq.(l3]), it is crucial to have the initial wave-functional (05’|in), 
equivalently the matrix element po[0i)02]- The goal of this paper is to derive the initial 
wave-functional for the perturbative Abelian gauge theory vacuum in the Milne space. An 
initial attempt in this direction was made in 1^. Within perturbation theory, non-Abelian 
gauge theory vacuum may be obtained by having A(? — 1 copies of the Abelian vacuum. We 
also verify that the propagators obtained in this way coincide with those obtained by solving 
the linearized equation of motion for the vacuum spectrum of fluctuations, the solutions 
being normalized with the help of the natural symplectic product QSi- 

In the context of relativistic heavy ion collisions, the initial state |in) is composed of the 
classical Yang-Mills held right off the light-cone and the zero-point vacuum huctuationsjl^- 


20(]. In this paper, we concentrate on the vacuum part. For the classical Yang-Mills part. 


see 


21 


25l |. The problem of vacuum huctuations in the classical background held is a more 
difficult one [l, 2| to which we hope to come back in a future work. 


II. KELDYSH ROTATION AND CLASSICAL-STATISTICAL 

APPROXIMATION 

The role of the vacuum density functional po[0i) ^ 2 ] can be better elucidated if one makes 
the Keldysh rotation 


0r — 


01 + 02 
2 


0a — 01 ~ 02 


(5) 

( 6 ) 
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Equivalently, 0i = 0r + 0a/2 and 02 = 0r~0a/2. The generating functional can be obtained 
by expanding £(0i) — / 1 ( 02 ) in powers of 0^ 


Ja\ 

= j V(j)r j P0aPw[0r,7r°] exp 


rTf 


drd^x ( (j)aE[(j)r] + 


(pl d^C{(j)r) 


24 


(/^0^ 4“ T^0f 


(7) 


where -E'[0''] = / 6C{(j)r)/S(j)r is the classical held equation, the canonical momentum held is 
dehned as tTj. = dC/d{dr(j)r), We also have performed the functional Wigner transform 

r r aO j .0 h 









( 8 ) 


where the 0a7rJ? term results from the boundary terms when integrating by part the terms 
containing 0r0a in E. The functional integration measure [^0°] denotes that the integration 
is performed with a hxed time, in this case at the initial time tq. 

The expression Eq.([7]) looks deceptively simple, but it needs a careful consideration before 
it can be useful. Notice that there is no integration over 7r° while pw depends on it. Recall 
that 7r° = dC{(j)r)/d{dr4>r)- Hence 7r° contains dr4>r evaluated at tq. This implies that pw 
is a functional of 0° and 0T-0r(ro). But how does one interpret the time derivative dr(j) in 
a path integral? This problem is not specihc to this particular formulation. In any path 
integral formulation of quantum mechanical system, the “time derivative” makes sense only 
in two ways. One is to discretize the time so that 0,-0 —)■ (0n+i — 0n)/Ar. The second way 
is to do perturbation theory using a classical solution as the zero-th order solution. 

When time is discretized pw[4>ry ^r] then provides restrictions on the hrst two held values 
0r(ro) and 0r(ro + At). This would be the case if one is to evaluate the path integral 
numerically. This is a fascinating subject by itself, but in this study we are more interested 
in the perturbative approach without discretizing the time. If we treat the 0j^£"Y24 term 
as a perturbation, then the unperturbed part of the generating functional is 


Zo[Jr,Ja\ = / [d0°][d7r°]p^y[0°,7r°] exp i / drd^X Jr4>cl[(pr,7^r,Ja] 


rTf 


' tq 


(9) 


using the fact that integration over 0^ produces (5[i?[0r] + Ja\ that enforces the classical 
equation of motion. Hence, 0ci[0Y7rY Ja\ here is the classical solution with the initial condi¬ 
tions given by 0° and vrj?. For more details such as converting the boundary value problem 


to the initial value problem, see Ref. 


26|. 
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If one can find an exact solution to the full nonlinear classical field equation, then the 
evaluation of an expectation value of any operator O can proceed by taking appropriate 
functional derivatives of Zo[Jr,</a] to match the expansion of in powers of 

C". In this perturbative approach, the time evolution is purely classical. Quantum effects 
are contained in the quantum fluctuations in pvi/[0°,7r°] and in the quantum correlation in 
the perturbing Lagrangian 0^£"'(0r)/24. This is the starting point of the classical-statistical 
approximation. 

If one cannot And an exact solution, then one needs to expand the interaction term in 
0a-E'[0r] as well. In that case, ZQ[Jr, Ja\ becomes the generating functional for the free-fleld 
and the resulting perturbation theory is the usual one with the free-fleld propagators. It is 
worth noting that the time evolution in a free-fleld theory is completely classical. Quantum 
effects reside solely in the fluctuations of the initial values. 

From Eq.Q, it is clear that the Wigner transform pw[(t>^r)'^r] can play the role of the 
probability density functional for the initial held values. Of course, as a Wigner functional, 
it is not guaranteed that pw is always positive definite. Fortunately as we show in this paper, 
the perturbative vacuum state usually results in a positive definite pw simply because it is 
usually a gaussian. It is also clear from Eq.Q that the two most important ingredients 
in obtaining Zo[Jr,</a] are the Wigner functional pw[(t>^r)^r] ^^cl the full classical solution 

(?!)cl[0°,7r°, Ja]. 


III. SUMMARY OF RESULTS 


Since this paper is of necessity technical, it will be useful to list the main results here 
before going into the details of calculations. Our gauge choice is Ar = 0. The spatial 
gauge vector field is further decomposed into the transverse vector and the longitudinal 
component A^. The 2-d vector A_l is further decomposed into AT{k) = k^A^ — k^A^ with 
k_i_ = k_i_//c_i_, and ip defined by V_l-A_i_ = 


The transverse component At satisfies the Klein-Gordon equation in the Milne space. 
Hence, we can adapt the results for the scalar theory found in Refs. 27|, l28|. The vacuum 
functional for At, re-derived in Section |Vl is known to be (omitting the 1,2 indices of the 
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Schwinger-Keldysh contour) 


{At \ vac) = N't exp 



d?k 

(27r)3 


AT{-~k) 



( 10 ) 


Here Nt is the normalization factor and (x) is the Hankel function of the hrst kind and 

of imaginary order, and we defined = xdxH[^’^\x). The Wigner transform of the 

density operator po = |vac)(vac| is in the retarded-advanced basis 




tTi '^rT. 


= Nt exp 


d?k 


RT{ToMKTm^ + 


|7r°y(fc) -F /r(ro, k)A%{k)\‘ 
Rt{to, k) 


( 11 ) 


where Rt and It are the real and the imaginary parts of —in(^(fc_i_ro)/-ffj^(fc_Lro), respec¬ 
tively. Eq. flTT]) is a new result. The Wigner function p^[A^rp,7i^rp] represents the zero-point 
motion of the quantum helds in vacuum whose effect manifests as the fluctuations of A^j, 
and which are the initial values of the ArT held evolution. It also plays the role of the 
probability density functional of A^j, and ir^rp in the context of the classical-statistical ap¬ 
proximation. In numerical implementations, Tq needs to be kept finite since neither Rt{tq, k) 
nor It^To, k) have a well defined limit as Tq —)■ 0. 

The r-a propagator is the retarded Green function of the d’Alembertian in the Milne 
space. It does not depend on p^. It is found in Section IVIIAI to be 


{ArT{r, -k)AaT{T', k)) = iGNr, k) 

The symmetric propagator is found to be 

{AMr, -k)AMr', k)) = | Nl{k^T)H^ll{k^T') + ff£>(fcxT)-ff™(fcxT')) (13) 

These propagators coincide with the Minkowski space retarded and the symmetric propaga¬ 
tors as shown in Appendix [Bl The exact same forms apply to the single component scalar 
held theory. 

For the longitudinal part, the two components and (p couples. The helds and 
ip themselves do not satisfy the free-held equation, but their conjugate momentum helds 
TT^ = (l/rjc^T-A^ and vr^ = —rV^dT-p do. Using this and the fact that the Gauss law requires 
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dr^T^n = it is shown in Section lYn that the vacnnm fnnctional is given by (omitting the 
1, 2 indices of the Schwinger-Keldysh contour) 


{Ari, (^1 vac) = Ml exp 




(14) 


where arj{k) = Arj{k) — ikjj(p{k). The Wigner transform of the vacuum density functional is 
found to be in the retarded-advanced basis 


Pw [-^r 


. tt: 


= Ml exp 


(Pk 

(^ 


IL(Ta,k)\AUk)\^ + 




rri V 


Miro, k) 


(15) 


where the real and the imaginary parts of k\H\^^{k±TQ )/have been denoted 
Rrii^o, k) and /,,(ro, k), respectively. Again, this can be interpreted as the probability density 
of initial values A^^{k) and 7r)?^(fc) which are the initial values for the held evolution. For 
the initial values of the (pr held, (p%k) can be set to zero using the residual gauge freedom 
and n^^{k) is given by ikr^7i^^{k). Note that having (pr{To, k) = pA^ik) = 0 does not imply (pi 
and p >2 both vanish at Tq. Only their sum is required to vanish. 

The longitudinal propagators in the r—a formalism are found in Section TVII Bl to be 

{Arnir, k)Aar,{T', k)) = -iOir - tq) f dr'A" (r'dr'CMr", r'; k) + toGr{t", tq; fc)5(r' - tq) j 

(16) 


and the symmetric propagator 


{Ar^{T, k)Arrj{T', -k)) 

^nS(fc±ro)ngi(fc±ro) ^ 7rfc2 


Akl 

r(2)| 

kr, ■ 


+ 


k)M^\T', -k) + ~k)M^\T', -k) 


-^(nS(’-o) [fl"(T,i) + fl"(T',-fc)] +n«(r„) [flM(ry) + flP'(T',-fc)]) (17) 


with 


R^'-^-HT,k) = / dTV Hl'P\k^T') 


( 18 ) 


'TO 


These are all new results. Now we turn to the calculation of these results. 
















IV. ABELIAN GAUGE THEORY IN THE MILNE SPACE 


In a curved space, the fundamental U{1) gauge field is a covariant vector Aa and the field 
strength is given by 


Fab = daAb - dbAa (19) 

The Milne space metric is given by 

ds^ = Qabdx^'dx^ = —dr^ + r'^drj^ + dx^ + d-ip ( 20 ) 

where t > 0, — 2 :^ > 0 while rj = z) can take any real value. The Abelian 

gauge theory Lagrangian in the Milne space is 

L = d^x C 

= _iy d^x^g^^g^^FabFad 

= - y d^XT {^rx + Ky + ~ ^xy ~ ~ (^1) 

where g = detgab = —and we set the coupling constant to unity. 

In the CPT formalism, the generating functional of the propagators is composed of the 
forward going time line and the backward going time line 


•^free ["-^l) "A] 

= /[dA-^][dA'j’][dA^]pvac[ro, A?, A^] [ VAi [ 

J Ja\ JaI 

( 22 ) 

Here Ai is the held on the forward going time line and A 2 is the held on the backward 
going time line. For obtaining the propagators only the quadratic part of the Lagrangian 
£free(A) is needed. The helds start with the initial values A\ and A® at the initial time tq 
and share the common hnal value A'f which is then traced. Note that the Lorentz indices 
and space-time dependence are suppressed in this expression for simplicity of notation. 

In terms of the ground-state wave-functional \ko[r, A°] = (A^l vac, r), the vacuum func¬ 
tional is given by 

= (A;| vac,r)(vac,r| Af) = Tolr, 


( 23 ) 
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Our task here is then to obtain \1/o[t, A'j’] that satishes the functional Schrodinger equation 


29 




idr^o[r,A]=H[A,T:]^o[r,A] (24) 

where H[A,ti\ is the functional Hamiltonian with the functional momentum operator tt = 
—i6/6A. Since the Hamiltonian turns out to be r dependent as shown below, this cannot 
be reduced to the time-independent Schrodinger equation. 

To obtain the functional Hamiltonian for an Abelian gauge theory, we hrst start with the 
Lagrangian in Eq.([2T]). In the A^- = 0 gauge, the Lagrangian density in Eq.([2T]) is (dropping 
the branch indices 1 and 2 for now) 

C = \ (VxX Ax)= + ^Ax-S^Ax 

+ 7^(^tA,j)^ -|- —- Arjdr^±-A.±_ (25) 

/r 2r T 

with the Gauss law constraint 

0 = + Vx-SxAx (26) 

We further decompose the transverse held as follows 


Aj_ = At + At (27) 

where 

^^^(k± X e^)AT(r, p, k^) (28) 

AT(r, x) = V_L(^(r, x) (29) 

These satisfy 


Vt-At = 0 (30) 

VtxAt = 0 (31) 

and 

/ (Px± At{t, x)-Al{t, x) = 0 (32) 

In terms of At and (p, the Lagrangian density becomes 

C = -(9.,-At)^ -I- -AtV^At -I- —AtO'^At 
l l It 

+ ’^{drArjY -l- —{dr'V+ — (A^V^A^j -|- V_L(p-5^V_L(p — 2Aj^cI,j(33) 
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The momentum fields are given by 


dC 


0^0 A X=^T = rdrAq 

0{drAT) 

^ = 7r„ = -d^A. 

d{drA^) 

dC 


r 


rj-J. J.^ 


^ TT^ ^ -rV\dr^ 


(34) 

(35) 

(36) 


and the Hamiltonian density is 


H = TTr^drAri + TT^dr^ + 'HtStAt — C 

1 7~ 1 

= — —AtV^At — —AtO'^At 

+ + 7^ (-A^VlA^ + (fd^Vlcp + 2A^d^Vlcp) (37) 


± 


2r 


The Gauss law constraint Eq. fl26|) becomes 




From the canonical quantization rule 

[Ha(r, (7), 7rb(r, ^)] = i5abS{x - y) 
we can identify the momentum-field operators 


ttt = —t 


At 0/1^ 5y) 

The functional Schrodinger equation then takes the following form 


vP- = -'W = ■* 


(38) 


(39) 


(40) 


idr'i>n[r,AT,A^,‘p] 
= H'i![T,AT,A^,ip] 


- I - ^ArVUr 

+ 7^ {— ArjW^Ar^ + + 2Aridri'V‘]_(p) j 'holG At, Arf, (p\ 


2 5AI 2t5^V\5^ 2t 


with the Gauss law constraint 




(41) 


(42) 


Note that since we are in the Schrodinger picture, the functional coordinates At, Ar,, (p does 
not depend on r. What we need in the generating functional Zfree[^i, ^ 2 ] is the solution of 
the Schrodinger equation evaluated at the initial time t = Tq. 
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V. TRANSVERSE PART OF THE VACUUM FUNCTIONAL 

Since the transverse and the longitudinal parts separate in Eq. lHTll . the wavefunctional 
factorizes 


d'oh At, An, tp] = ^t[t, At]^lIt, An, (f] (43) 

We consider the transverse part hrst in this section. The longitudinal part will be discussed 
in the next section. 

The time dependent Schrodinger equation Eq. fl42|) is valid for any quantum state, not just 
the vacuum state. If the Hamiltonian is time independent (for instance in the Minkowski 
time t), then we can let 

A] = e [A] (44) 

and solve for the time independent equation 


E4/t[A] — i^Mink.[A, 7r]\k'r[A] 


(45) 


However in the Milne space, evolution in r does not preserve total energy since the Hamil¬ 
tonian in Eq. dTT]) is explicitly r-dependent. Therefore, the usual concept of vacuum as the 
lowest energy eigenstate of the Hamiltonian does not apply. 

On the other hand, the Milne space is just a re-parametrization of the flat Minkowski 
space. Therefore, the state representing vacuum should have a similar structure as the 
Minkowski vacuum functional. In the Minkowski space, the perturbative vacuum for a 
scalar mode is known to be 


/ If ^ 3 ^ \ 

d'Mink.lt, A] = exp - 2 J 

where Eq = V j E^d^k/{27r)^ is the un-renormalized vacuum energy. The corresponding 
state in the Milne space should also have a Gaussian structure. Furthermore, for the two 
space-time points inside the Milne space, they should result in the same two point correlation 
functions (propagators) as in the Minkowski space. Hence, we employ the following Gaussian 
ansatz for the vacuum functional 


^t[t. At] = Mt{t) exp 


2 


d?x 


d^y At{x)Gt{t, X 


y)AT{y) 


(47) 
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Since we are in the Schrodinger picture, At{x) does not depend on r but the kernel Gt{t,x — 
y) does depend on r. Consistency under x y exchange requires 

Gt{t,x- y) = GT{T,y - x) (48) 

By comparing the left hand side and the right hand side of the functional Schrodinger 
equation (Hill , one obtains the equation for the normalization factor 

i9^A/'r(r) = ^ J d^xGT{T,0) XTir) (49) 

and the equation for the transverse kernel 

idr{T‘^GT{T, X — y)) = J (fuGrir, x — u)Gt{t^ u — y) + H— 5{x — ^)(50) 

The normalization factor obviously involves inhnities since the integrand does not 

depend on x. It requires renormalization. However since Mrix) is independent of the field 
variables, it does not concern us in this paper. 

In the momentum space, the equation for Gt simplihes 

~ / k‘^\ 

idrir’^GTir, k)) = t^Gt{t, k)GT{T, k) - i Tk]_ + j (51) 


Following the Refs.j^. 28|, dehne 

- idri>Tir, k) = tGt{t, k)‘ipT{r, k) 

Multiplying 'ipT{,T,k) to Eq. fl5T]) and using Eq. flS^ . Eq. fIFT]) becomes 


(52) 


= 0 


(53) 


where 


= a! 


-a. 

r 


+ k] 


N 

r2 


(54) 


is the d’Alembertian in the Milne space. 

The solutions of Eq. flS^ are the Hankel functions of imaginary order if^^(fcxr) and 
Hll^^{k±T). If vacuum is to be stable, the real part of Gt must be positive dehnite. Further¬ 
more, it must be able reproduce the Minkowsi correlator when the two points are both inside 
the forward light cone. For these to hold, we must choose iprix, k) = cH\^^{k±_T) where c is 
a constant. Using Eq. (|5^ one obtains 


t^Gt{t, k) 




(55) 
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The real part of k) is indeed positive definite. 

Rt{t, k) = Re(r^GT(r, k)) 

= -rp-^-> 0 (56) 

where we used the Wronskian of the Bessel functions 

Hlt’(.T)drHS\T) - H‘i(T)drH[^(T) = (67) 

TIT 

and the fact that H‘^^{k±T) = e'"^"'{Hll^^{k±T))* (see Appendix IXI for Hankel function prop- 

/'o\ ~ 

erties). Had we chosen H-f^^{k±T) as our solution, the real part of Gt{t, k) would have been 
negative definite. It is convenient to define 

nSif (fc±r) = TdrH\lf{k^T) (58) 

as they occur frequently in our analysis. The fact that ^^[t, A^] is compatible with 
d^Mink. A't] is demonstrated in Section IVlI Al bv showing that the resulting propagators 
are the same. 

The matrix element of the vacuum density operator at tq is now given by 


Pvac [ A)) A 


0 

IT) 



= T' 


^[to, A^’Ad' 


= exp- 


T [A) ) ^ 

d?k 

W) 


° 1 

2T\ 


A^rp{ 


-k)GT{To, k)Al^{k) + AlT{-~k)G*T{To, k)AlT{k) 


(59) 


VI. LONGITUDINAL PART OF THE VACUUM FUNCTIONAL 


For the longitudinal part, our ansatz is 


T 


TiiA A^,(p] = A/'L(r)exp I -y / d^x / d^y [Ar^{x)Gr,{T,x - y)A^{y) + (p{x)G^{T,x - y)(p{y) 


+^p{x)G^n{T, X - y)Arj{y) A^(x)G^^(r, x - y)^p{y)] 


( 60 ) 
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Since the normalization factor Ml does not depend on the fields, again we will not consider 
it. Consistency under x ^ y exchange requires 


G,(t, X - y) = G^(T,y - x) 
G^{t, X - y) = G^(T,y - x) 
G^(t, x-y) = G„{t, y-x) 


(61) 

(62) 

(63) 


Applying the Gauss condition Eq. (l42p on the ansatz, one obtains the following relationships 
between the longitudinal kernels 

d 


dr], 

d 

dr]x 


x-y) = x - y) = y-x) 

X 

G„(r,i - y) = G^(t, x- y) 


(64) 

(65) 


Comparison between the left hand side and the right hand side of the functional Schrodinger 
equation yields the equation for the longitudinal kernel 

- «5.(r^G,(r. x - y)) ^ J SaG,(r, x - i)G,(r. u - y) 

j3 ~ ^3- 


+ r'^ y d^u (TV x — u)Gl{u — v — y) 

- -Vi(5(x - y) 

T 


( 66 ) 


where 


G_l{u-x) = S{r]u-r]x) 


V2 


(67) 


The equations for G^ and Gn^p are the same as Eq. (l68l) upon applying the Gauss constraints, 
Eqs. dM]) and fl6^ . In the momentum space, the equation for becomes 

1 


where 


k‘^\ 

idX,,{x, k) ^ T [ k'i + ^ j k)0^(r, ^ 


Sy(T,k) = ^G,,(T,k) 


( 68 ) 


(69) 


and we used Eq. (|64|l to express G^,, in terms of G^. Let 




- idMr^ir, fc) = r ( fc_L + ^ ) i/’^(r, k)Qnir, k) 


( 70 ) 
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-dr{TdrX{r, k)))+ [k\ + ^] x(r, fc) = 0 


(72) 


and multiply Eq. fl68|) with to get 

- irdriGnir, k)'ijj^{T, k)) = k) (71) 

Combining Eg. (1701) and Eg. (17X1) . one gets 

1 

— ( 

T 

where xij^k) = 'ip{T,k)Qr]{T,k). This is again the free held equation for y. Stability of 
vacuum again selects Hllj^{k±T) 

~k)G^{,T,k) = c'H\^^{k^T) 

where d is a constant. Putting this back into Eq. fITT]) gives 

i\)n{T~k)/d = -iTdriif^^{k^T) = -ilif^^ik^T) 

so that 


(73) 


(74) 




T‘^Gn{T,k) = 




(75) 


Our choice of the solution ensures that the real part of the kernel is positive 

dehnite 


7?^(r, fc) = Re (r^G^(r, k) 


''1. 
2 

2fc2 


r(i) 

kr) 

r(b 




^Zik±r) IlZik±T) 


( 2 ), 

7] ' 


> 0 


^ n«(fc^T)nP|;(fcj_T) 

The longitudinal part of the vacuum density functional is now 

Pl [t'O) ^ini 7^15 442j,, P2y ] 

= (7l?,,v'!|vac,ro)(vac,ro[/l2,,v>2) 


(76) 


= M exp — 


Tn 


d^k 

(^ 


a°v(-pG^(Ta, k)al^(k) + a5,(-(:)G*(Tci, k)al,^(k) 


(77) 


where we used the Gauss constraints Eqs. fl64l) and fl65ll to express it solely in terms of the 
combination a)* 2 ,, = ^\, 2 ri ~ ^^e momentum space, 

“i,2,(^) = X,2„(k) - ik^ipl^iX 


(78) 
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VII. WIGNER TRANSFORMS AND PROPAGATORS 


A. Transverse part 

With the vacuum wavefunctionals at hand, we can now calculate the propagators. The 
calculation is simpler if one changes variables to the r-a basis (also called the Keldysh 
rotation) 


ArT — {Ait + A2t)I‘^ 
AaT = Ait — ^2r 


(79) 


because the separation of the classical and the quantum degrees of freedom is cleaner in this 
basis and also because the AaT—AaT correlator automatically vanishes. 

In this basis, the transverse part of the generating functional Eq. fl2^ becomes 


t\ [^^ar] Pvaclh)) A^j,, 



VA 


rT 


VA 


rT 


'Kt 




X exp [i dr d?x ( AaT 


SCxiArT) 


5A. 


+ JjArT + jjAaT ) “ M d^^TT^TAaT 


rT 


where A^j^ = Arxiro), A^j, = AaT{To) and 

T^rT = T^rT{To) = 


dCfree{ArT) 


(80) 


(81) 


•^0 


d{drArT) 

is the momentum held at r = tq. The term containing in Eq. flHOj) came from integrating 
by part in r. As in Eq. fl2^ . we have suppressed the space-time indices above. 


^•^free (-^rx) 


6ArT 




which enforces 


Performing the j VAaT integrals for r > tq results in <5 
the classical held equation. Performing the integral at tq produces the functional 

Wigner transform of the vacuum matrix element 

d^k 


PwivArT^T^rT] = / [dA^-^] Pvac[E ^an] eXp -f 


(2vr)- 


AU-kKTik) (82) 


It can be shown that the Jacobians coming from integrating over the d-functional and chang¬ 


ing integration variable from A{.rp to tt^t (transformation between the bounc 
lem and the initial value problem in classical held theory) cancel each other 


ary value prob- 


261. Hence, the 
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generating functional becomes 

ZtIJ't, Jt\ 

^ exp[i j dr / J^A^rl^rr* ^rTi 


(83) 

The transverse field ArTiA^lj., J^] is the solution of the equation of motion 

Td'^ArT + Jt = 0 (84) 

with the initial conditions ArT^To) = ^4°^, UrTiTo) = TodrArriTo) = 7r°^. In the momentum 
space, the full solution is given by 

poo 

Arrir, k) = A^t{t, k) - I dr Gr{t, r'; k) k) 


where Gr is the retarded Green function defined by 




The solution is 


Gr{t,t’, k) = -0 {t - t') - ffS>xr)i/S>(A.ir') 

which can be easily verified by the direct substitution into the defining equation 
using Eq. fl57|) . Performing Fourier transformation in kr^, one obtains (see Refs, 
also Appendix iBl) 


Gs(r, k±) = -0(7 - T’)e{a^)Jo(kRV^) 


(85) 

( 86 ) 

(87) 

and 
33j and 

( 88 ) 


with cr^ = (t — — {z — z'Y = — 2rr'cosh(?7 — rj') > 0. This coincides with the 

Minkowski space retarded propagator. The r-a propagator can then be read off from the 
J^—J^ part of the exponent in Eq.(jH3]) as 

{Arrir, k)AaT{r\ x = 0)) = iGRij, r', k) (89) 

To obtain the symmetric propagator {Arrir, k)ArT{T', —k)), one needs to explicitly evalu¬ 
ate Eq. dH^ . To do so, we need the homogeneous solution A^j,{t, k) and the Wigner transform 
Pw[ro, With the retarded Green function given in Eq. fIST)) . the source term van¬ 

ishes in the r —)■ tq limit for both and drA^T- The homogeneous solution satisfying the 
right boundary conditions is 


(90) 


+ n (n® (A^xT„)ff<,i;(A^xT) - n™(A^xT„)ffg;(A^xT)) A%[k) 
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Using the Wronskian Eq. flHTll . it is easy to check that ArxiT, k) given in Eq. fl^ does satisfy 
the initial conditions. Since the held in the conhguration space must be real, ~ 

{T^%{—k))* and A^j,{k) = [A^rp[—k))*. It is trivial to check that the above solution also sat- 
ishes A^j,{t, —k) = {A^j,{t, k))* using the Hankel function property 

In the r-a basis, the vacuum density functional pvac[To, ^ 2 t] Eq. (j5^ becomes 


Pvac[ro, A%, AIt] = exp - 


1 

2 


(Pk 

{2nf 
'• d^k 

W) 


Rx{To,k) 


■ kriro, k) 


AP^t{ k)A^^rp{k) + .^aT( ^)^ar(^) 




-~k)A%m + A^,^{-k)A%{k) 91) 


where RxijQ^k) = Re(rQ(7^(10, fc)) and IxijQ^k) = ImirQGxijQ^k)) with x^GxirQ^k) = 
—in|^(/c_i_ro)/iff^(/c_Lro) as given in Eq. flS5]l . The functional Wigner transform can be 
carried out by completing the square in the exponent 

f d?k 


PwVo^Kt^t^tx] = / [dKr] 


(27r)- 




= M exp 


d^k 


Rx{To,k) A%{k) 


+ 


7r°r(^) + lT{rQ,k)Alj,{k) 


Rt{to, k) 


(92) 


where M is the normalization constant. This can be interpreted as the probability distribu¬ 
tion of the vacuum huctuations of the helds A^j, and at tq. 

The generating functional for the symmetric propagator can be evaluated by completing 
another square 


Zs[Jx] = I [dA^^j][d7i^T]pw[To,A^^j,,Ti^j,] exp ( i / dr I 


'TO 


= AA'exp -- 


drdr' 


d^k 


J^(r,-fc)J^(r',fc)Gs(r,r';fc) 


(93) 


where J\f' is another normalization constant. The symmetric propagator is 
{Arxir, k)Arx{,T\ X = 0)) = Gs{r, r'; k) 

+ TfS^(mrr)i/£J(mTr')) (94) 


TT 


Remarkably, the dependence on To has gone away. Similar result was obtained in 3^. Fourier 
transforming to (r, 77,kx) space yields (see Appendix 1X1 for Hankel function properties) 

Gs(t,t',‘I]- ri',kx) = -S{a^)^Yo(mT\XX) + 9(-<T^)XA'o(m7-\/-ff2) 


(95) 
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This is the same as the Minkowski space symmetric propagator, further justifying our choice 
of the solution for 'i/lr(r, fc) = in Section |Vl The same vacuum spectrum and 

propagators apply to the single component scalar held theory. 


B. Longitudinal part 


After the Keldysh rotation 


Oirr/ T Cl2r;)/2 

O-ar) ^Ir} ^2ri 

the Longitudinal part of the vacuum density functional Eg. (17711 becomes 


(96) 


Pl[to, = Af exp 


d^k 


(2vr) 


Rn{To,k) 


al^{-k)al^{k) + jal^{-k)al^{k) 


d^k 


4(^-0, k) 


a 


ar)\ 


k)a:° (k) + a° (-k)al (k) 


(97) 


where k) = Re(rQG'^(ro, k)) and /(^(ro, k) = Im(rgG^(ro, k)) with TqG^ as given in 

Eq.(I75D. 

To perform the functional Wigner transform, we hrst add the boundary terms to the 
exponent 


d?k 


— i 




d^k 


«ar,(-^)^rr,(^) (98) 


(27r)3 V rrjX-'J ' TTarjV V " J (^27ri^ 

d the Gauss law constraint and 

the square then yields the Wigner transform just as in Eq. 


where we used the Gauss law constraint and integrated by part. Gompleting 


Pw[^o,a™,7r° ] =U' exp 


d^k 


Rn{To,k) (k) 


+ 


7r° (fc) + 4(ro,fc)a° (fc) 


2"l 


Rniro^k) 


(99) 


This has the interpretation as the probability distribution of the vacuum huctuations of the 
helds and at tq. 

Since depends only on the combination — ikrj^p^ individual propagators 

cannot be calculated yet. Fortunately, there is still one more residual gauge freedom left in 
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the A,- = 0 gauge. That is, the following gauge transformation is still allowed 

= A±{t,x)+ V±f{x) (100) 

x) = Ar,{T, x) + dr,f{x) ( 101 ) 

as long as f{x) is independent of r. Since V_l-A_i_ = the effect of Eq. fllOOl) on tp is 

ip'{T,x) = p{t,x) A f{x) (102) 

The momentum held vr^ is unaffected by this transformation. We can choose f{x) = —pP{x) 
so that 

(p'{to,x) = 0 

and 

^'r,(T'o,5) = a7,(5) 

In this way, depends solely on and 7r°^. 

To obtain the longitudinal propagators, we start with the generating functional containing 
the vacuum density functional 

J^] = j [d7r°^] p^[ro, vr)?^] exp J J (105) 

where we now dropped the prime on A^rj- To obtain the propagators, one needs to hnd the 
solutions of the classical held equations 

rdr (^drAr^ - V]_Arrj + = 0 (106) 

Tdr {rdrWlpr) " + Vld^^Ar^ - T= 0 (107) 

with the initial conditions determined by A^^ and Although A^r^ and pr themselves 
cannot be solutions of the free-held equation, the momentum helds dehned by 

Til - ^T-^^TTj 

T 

T^np = —rdrV\pr (108) 

do satisfy the free-held equation 

/I d‘^\ 

T ( + -dr TTr,,,^ + dr(Tj“J = 0 (109) 


(103) 

(104) 
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upon using the Gauss law constraint d^Tirri = The sources must also satisfy 


8 r = 7°- 


(110) 


again due to the same Gauss law constraint. 

The solution for Tr^rj is given by 

T^rrjir, ^ 7r,^(ro,fc) 

+ ^ + Hlll{k^To)Hj;^l{k^T)'^ (roa,7r,^(ro, k)) 


JTQ 

ivr 

4 V ‘"■’7 
ITT 


dr' Gr{t, r'; A;)9^(r'J“(r', k)) 

Il[ll{k±To)Hll\k±T) - u[lHk^To)Hl^\k±T)) 7r,^(ro, k) 


+ -J- [--^lk!,(k±'^o)kf!l!,(k±T) + Illl^(k^To)I7l^^(k±T)j (roa,7r,^(ro, fc)) 

+ / dr'(r'9^/GH(r, fc)+ roGH(r, ro;fc)J“(ro,fc) (111) 


where the retarded Green function Gr is given in Eq. flHTl) . For TTr^, the solution is simply 
k^ xk'qTly^^i^'T^ A)). 

The helds are given by 


A()(r, k)= dr' r' 7rr^(r', k) + Al^{k) 


( 112 ) 


' ro 


and 


( r\ r dr' ik,j7lrrj{T',k) 

^r{r,k) = ' 

J TO 




(113) 


Putting this back into the equations of motion produces the following relationship between 
and dr'llrr, 


= -j2TATrn(To,k) 


(114) 


As before the r-a propagator can be read off from the classical solution Eq. flll2p as 


(A^^(r, k)Aar,{T', X = 0)) = - 


6 


-A^jj(t, /c) 


idJair', k) 

= -i9{T - To) / dr"r" ( r'dr^GR^T", r'; k) + toGr{t”, Tq; A:)5(r' - Tq) 

(115) 


' ro 
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where from Eq. fl87p 


TT 


r'5.,G„(r",T';«:) = ^«(r" - r') (4‘>xr")nS'(A^ir')) (116) 

The r-a propagator for the (p field can be similarly obtained using 

dr" 




r(2) / 




((p^(r, k)ipa{T', x = 0)) = - ^o) 


'TO 


r" 


T'dr'CR^T", r'; k) + roG'/j(r", tq; fc)h(r' - tq) 


(117) 


To obtain the symmetric correlator, we need the homogeneous part of Arrj from Eq. flll2p 




J“=0 


= K(r, + T(r, k)7r^^{k) 


r?7 V 


(118) 


where 


and 


with 


/^(r, fc) = 1 + 


ZTr/c^ 


-L ( tt {‘ 2 ) , 


HZ‘.kLT,>)R?{T, k) - ii;;'(l;xTo)-Rr'(G k) 


'b 


,(2), 


i(T, ^ “ n,(s!,(*^i^o)-RfT. k) 


(119) 


( 120 ) 


/ <lTV“J/Of(4xr') 

J TQ 

Putting this into Eg. (11051) and completing the square, one gets 


{Arr^{T, k)Ar^{T',X = 0 )) 

^n2(/c±ro)n® (fc^To) 7rfc2 


4A;2 


+ ^ ( R? (r, k)R\^> (r', k) + R\^> (r, k)R\^> (r', fc) 




( 2 ), 




^ (n(^)(ro) [R^^\T,k) + R?{T',k)\ +n«(ro) [Rf\T, k) + R^^>{r', k)\) (122) 
The TTrrj—'^rrj Correlator can be obtained by applying {l/TT')drdr' to {Arrj{T, k)Arrj{T', —k)): 
(7r,^(r, fc)7r,^(r',x = 0)) = ^ fi7£VA:xr)i7£(A)±r') + id® (/cxr)i7£(/c^r')) (123) 




( 121 ) 


The (fr—^r correlator can be obtained using Eq. flll3p 


nk‘? 


{cprir, x)ipr{T', X = 0)) = ^ [ R^}1{t, k)R^^K t , k) + i?'Lj’(r,/c)i2!_'((r',/c)) (124) 


(2)^ / 


>(2), 
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VIII. COMPARISON WITH THE SYMPLECTIC PRODUCT METHOD 


Another method to obtain the propagators amounts to solve the linearized equation of 
motion in order to obtain the vacuum spectrum of fluctuation^ the solutions then being 
normalized with the help of the natural symplectic product jl, 2|. To see if the symplectic 
product method is compatible with our approach, we re-arrange Eq. fMI]) as 


ITT 


KtM = y ( -HSy±nKT{>=) + ^'£,{k±n)A"^{k) ) 




't’l 


+ J - n‘il{k^T,)AUk)) 


(126) 


The vacuum functional in Eq.([92]) has the exponent 


St = - 


d^k 

(27r)3 


A I^T A /Rt) 


(126) 


where Rt and It are the real and the imaginary part of the kernel TqGt{to, k) as given in fl55l) . 
Since A^rp{x) and are real, we have A^j,[—k) = [A^rpik))* and 7rJ?y(—fc) = (vr° 2 .(fc))*. 

Hence the real part of AP^rp{k) and Ti%{k) must be an even function of k and the imaginary 
part must be an odd function of k. The real and imaginary parts of TqGt{to, k) are also 
even functions of k because H^^f,^{k±T) = (see appendix section |^. Hence 

under the J d^k integral, we can replace 


Rt\A 


^ ^ \7r'^j,-i{RT + ilT)A^T\ 

rT\ “I ^ 


Rt 


Rt 


(127) 


because the difference is an odd function of k. This allows us to re-express the exponent in 


w 

prp as 


St = -2 


^ ^ '7r°;r(fc)ai(ro) - A°j,{k)el{To) 


(27r 


(128) 


where used the fact that 1/Rt = '^Hil\k±To)Hl‘^\k±To)/2 and defined 




(129) 


and 


ej(T) = Tdra^(T) 


(130) 
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For large r, a^(T) 






\/2k±T 

Hence a^(r) can be regarded as the positive frequency solution. Since 

the Wronskian Eq. flSTl) is equivalent to the symplectic product 

e*_k{T)a-^{r) - e-^{T)a*_-^{T) = i 
In terms of aj^ir), the homogeneous solution is 

^rT(F k) = c^a-^ir) + 

where 

Cfc = i (^%(k)a*k^{k^To) - H°r(fc)e.^(/cxro) 

In terms of c^, the exponent of becomes much simpler 

d^k 


St = -2 


(27r 


,3 ‘^k'^k 


which implies 


(Cfc4) = 


0 


lis is identical to the correlation function obtained with the symplectic product 
2|- 


For the longitudinal part, the exponent in is 


S,= - I SS(Rr,\A«„P + h'” ^ 


Rn 


J (2vr) ^ 

where and is the real and the imaginary part of ik‘j_Hll^^{k±TQ){k±T). 
be re-expressed as 


Sl = 


d?k |vr°^ - 


"n/ r?jl 


(27r) 


R-n 


(131) 

(132) 

(133) 

(134) 

(135) 

(136) 

(137) 
approach 

(138) 
This can 


(139) 
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again using the fact that the real parts of and are even functions of k and the 
imaginary parts are odd functions of k. In terms of afc('r) and e^(r) dehned in Eqs. fll29p 
and (11301) . we have 

2 


Rr, + Hr, = ik 


ey(ro) 


and 


Hence 


Sr = 


k‘2 

2efc(ro)e;^^(ro) 


dPk \7r^ —i{R,j +Hr,) A^. 


TT] I 


= -2 


(27r)3 

r (Rk 1 


Rn 




2 _ /lO 

k^rr}\ 


j {2 tt)^ k\ 

The same combination in the integrand occnrs in (c.f. Kn. lll i m i 

tre,(T, k) = diaj(T) + d* ^a* ^(r) 


(140) 


(141) 


(142) 


(143) 


where 


d'k = R-AA(k) + k\al{To)AA) (144) 

Again, the exponent simplihes in terms of R 

/ Rk 1 

which results in 

(d-kRy) = Y^ARk - k') (146) 

which is again equivalent to the correlation function used in the symplectic method. 

Two remarks are in order here. The hrst one concerns the number of independent stochas¬ 
tic degrees of freedom. If one compares Eq. fll26p and Eq. (11361) . it looks like the number of 
stochastic degrees of freedom is reduced by a factor of 2 between Eq. (ll26p and Eq. (ll36p . This 
is not so. Even though it looks like there are 4 functional degrees of freedom in Eq. (ll26p . 
they are not complete. Recall that the real part of A^^rp and is an even function of k 
and the imaginary part is an odd function. Thus, each cover only 1/2 of the function space. 
On the other hand, there is no such restriction on in Eq. (ll35p and hence there is no 
contradiction. The second remark is that for numerical implementation Eq. (ll36p is much 
more suited than Eq. (ll26p . which is what has been used in 351. 
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IX. CONCLUSION 


In this paper, we have derived in Milne space the initial density matrix of an abelian gauge 
field theory in vacuum, by solving the corresponding Schrodinger equation. This allowed us 
to compute the Wigner transform of the density matrix, leading to a simple expression for 
the transverse and longitudinal parts of the generating functional of the Closed Time Path 
formalism. From there can finally be deduced the propagators of the theory. The key finding 
of this study is the perfect agreement between the results obtained with the functional 
Schrodinger approach and those derived with the symplectic product based method - 
as demonstrated in section IVIIII An important though challenging extension of this work 
would be to compute the initial density matrix of the Wave-functional corresponding to the 


so-called Glasma state 


21 


23| . and to compare with the result of in this interacting case. 


Establishing with the Schrodinger functional approach the Glasma initial condition is 
of crucial importance for dealing with the problem of thermalization in heavy ion collision. 


IndeeT 


fluid 


36 


le quark-gluon plasma created in these collisions seem to behave like a nearly perfect 


37| on very short time scales, as corroborated by numerous hydrodynamical studies 


that are able to reproduce the experimental data collected at the RHIG and LHG colliders. 


But a hrst principle proof of the app 
- that is initially not thermal}?. 


icability of hydrodynamics for the quark-gluon plasma 
I - has so far been lacking. An important building- 


block towards the completion of this proof is the interacting counterpart of the present 
work, derived with the symplectic product based approach j^, that has recently been used 


numerically to obtain promising results 


Improving our theoretical understanding of 


the Glasma initial condition with the Schrodinger functional approach would therefore be 
of great help in the path towards the full understanding of the nearly perfect fluidity of the 
quark-gluon plasma. 

Another challenging extension is the NLO contribution to the vacuum wavefunctionals. 
At the leading order, there is no real difference between the Abelian vacuum and the non- 


Abelian vacuum. However, the NLO contribution must be differen 


developments on the NLO JIMWLK equations and factorizations 4fll-l44|. this would be an 


In view of recent 


important step to establish firm connection between the current Schwinger-Keldysh approach 
and the GGG approaches. 
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Appendix A: Hankel function properties 


The dehnition of the Hankel functions used here are 


A‘A) = 


• 7rfc,,/2 poo 


le 


n 


dC e 


iX coshCT^GC 


(Al) 


and 


r) ^ 


— l2 rOO 


dC coshC±ikr}C 


(A2) 


with real y. Both signs of ikr/C leads to the same function because they are related by the 
simple change of the integration variable from ( to From these dehnitions, it is easy to 


see 






-,-TTkr, TjW 


( 2 ), 


(A3) 

(A4) 


Using the dehnitions, one can perform Fourier transformation of the product of Hankel 
functions. Let 

Inir, t', 7]-7]')= f ^ 


= — / dCe 


imr cosh((^— 77 ) —imr' cosh((^— 77 ') 


TT^ 


' —CX) 

poo 


= ^ dCe 


cosh C, — {z—z') sinh^) 


TT^ 


(AS) 


where t = rcosh?] and = rsinh?]. If [t — t'Y > {z — z'Y, then the exponent can be 
combined into 


{t — t') cosh — {z — z') sinh C = sign(t — t')a cosh(C — Co) 


(A6) 



where Co = tanh ^{{z — z')/{t — t')) and 


a = \/ {t — t'Y — {z — z'Y 

= a/t^ + r '2 — 2 rr' cosh (?7 — 77 ') 


Hence 


luir, t',7] 


Similarly one can show 


r';v-v') 




j ^ gifc .(V V)^(^2) ^ w 


and since Hq^\x) = Jo{x) + iYo{x), we get 

, 2 

Ih + Ih = - yo{ma) 

71 

as well as 

2i 

Ih-I'h = 

71 


(A7) 


(AS) 


(A9) 


(AlO) 


(All) 


Appendix B: Transformation of Minkowski space propagators 

We want to show the equivalence between the symmetric and retarded propagators in 
Milne space and the usual Minkowskian ones. Starting from the symmetric propagator in 
the Minkowskian coordinate system 

Gs{k) = TTdikl - El) (Bl) 

we can re-express it in the following way, depending on the region of the light cone considered 

(1) In the forward light cone region, where t^—z'^ > 0 and t > 0, a straightforward calculation 
leads to 

Gs = -^lo(mTr) (B 2 ) 

(2) In the backward light cone region, where — 2 :^ > 0 and f < 0, one easily hnds the same 
result as in the forward light cone region 

Gs = -^lo(mTr) 


(B3) 
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(3) In the right side of the light cone, where and z > 0, one can use the fact that 


to show that 


k^z ± E^t = mTcr sinh.{y ± 7]) 


Gs = —KoimTT) 

/TT 


(4) In the left side of the light cone, where and z < 0, one also hnds 

Gs = —Ko{mTT) 
zn 

which gives back fl95ll as anticipated. 

Now for the retarded propagator, we start with 

sm{Et) 


GR = e{t)- 


E 


which can be re-expressed into 

Gr = 

Consider the following two cases 
(1) If we assume \z\ > |f|, then 


dk^ g— z—iEt 

2iE 


k^z — Et = sign(z)m'r\/—cr^ sinh(|/ — sign( 2 ;)r 7 ) 


leads to 


/ dy 

— sin(a;sinh( 2 /)) = 0 

(2) If we now assume |f| > \z\, then 

k^z — Et = —mTV^sign(f) cosh(?/ — sign(t)r]) 

leads to 

Gr = 


(B4) 


(B5) 


(B6) 


(B7) 


(B8) 


(B9) 


(BIO) 


(Bll) 


(B12) 


as announced in 


[1] K. Dusling, F. Gelis, and R. Venugopalan, Nucl.Phys. A872, 161 (2011), 
arXiv:1106.3927 [nucl-th], 







30 


[2] T. Epelbaum and F. Gelis, Phys.Rev. D88, 085015 (2013), arXiv:1307.1765, 

[3] E. lancu, A. Leonidov, and L. McLerran, “The Color glass condensate: An Introduction,” 
(2002), 73-145 pp., arXiv:hep-ph/0202270 [hep-ph], 

[4] A. Kovner, Acta Phys.Polon. B36, 3551 (2005), arXiv:hep-ph/0508232 [hep-ph], 

[5] J. Jalilian-Marian and Y. V. Kovchegov, Prog.Part.Nucl.Phys. 56, 104 (2006), 

arXiv:hep-ph/0505052 [hep-ph], 

[6] H. Weigert, Prog.Part.Nucl.Phys. 55, 461 (2005), arXiv:hep-ph/0501087 [hep-ph], 

[7] T. Lappi and L. McLerran, Nucl.Phys. A772, 200 (2006), arXiv:hep-ph/0602189 [hep-ph], 

[8] L. McLerran, Acta Phys.Polon. B41, 2799 (2010), arXiv: 1011.3203 [hep-ph], 

[9] F. Gelis, E. lancu, J. Jalilian-Marian, and R. Venugopalan, 

Ann.Rev.Nucl.Part.Sci. 60, 463 (2010), arXiv: 1002.0333 [hep-ph], 

[10] 1. Balitsky, “High-energy amplitudes in the next-to-leading order,” (2010), 
arXiv: 1004.0057 [hep-ph], 

[11] Y. V. Kovchegov, AIP Conf.Proc. 1520, 3 (2012), 

[12] E. lancu, “QCD in heavy ion collisions,” (2012), arXiv: 1205.0579 [hep-ph], 

[13] F. Gelis, Int.J.Mod.Phys. A28, 1330001 (2013), arXiv:1211.3327 [hep-ph], 

[14] J. S. Schwinger, J.Math.Phys. 2, 407 (1961). 

[15] L. Keldysh, Zh.Eksp.Teor.Fiz. 47, 1515 (1964). 

[16] K. Fukushima, F. Gelis, and L. McLerran, Nucl.Phys. A786, 107 (2007), 
arXiv:hep-ph/0610416 [hep-ph], 

[17] F. Gelis and R. Venugopalan, 

[18] F. Gelis and R. Venugopalan, 

[19] F. Gelis, S. Jeon, 
arXiv:0706.3775 [hep-ph], 

[20] F. Gelis, T. Lappi, and R. Venugopalan, Phys.Rev. D79, 094017 (2009), 
arXiv:0810.4829 [hep-ph], 

[21] L. D. McLerran and R. Venugopalan, Phys.Rev. D49, 2233 (1994), 

arXiv:hep-ph/9309289 [hep-ph], 

[22] L. D. McLerran and R. Venugopalan, Phys.Rev. D49, 3352 (1994), 

arXiv:hep-ph/9311205 [hep-ph], 

[23] L. D. McLerran and R. Venugopalan, Phys.Rev. D50, 2225 (1994), 


Nucl.Phys. A779, 177 (2006), 

arXiv:hep-ph/0605246 [hep-ph] 


Nucl.Phys. A776, 135 (2006), 

arXiv:hep-ph/0601209 [hep-ph] 


and R. Venugopalan, Nucl.Phys. A817, 61 (2009), 





31 


arXiv;hep-ph/9402335 [hep-ph], 

[24] A. Ayala, J. Jalilian-Marian, L. D. McLerran, and R. Venugopalan, 

Phys.Rev. D52, 2935 (1995), arXiv:hep-ph/9501324 [hep-ph], 

[25] A. Ayala, J. Jalilian-Marian, L. D. McLerran, and R. Venugopalan, 

Phys.Rev. D53, 458 (1996), arXiv:hep-ph/9508302 [hep-ph], 

[26] S. Jeon, Phys.Rev. C72, 014907 (2005), arXiv:hep-ph/0412121 [hep-ph], 

[27] D. Long and G. Shore, Nucl.Phys. B530, 247 (1998), arXiv;hep-th/9605004 [hep-th], 

[28] D. Long and G. Shore, Nucl.Phys. B530, 279 (1998), arXiv:gr-qc/9607032 [gr-qc], 

[29] J. Greensite, Nucl.Phys. B158, 469 (1979), 

[30] A. Jevicki and B. Sakita, Nucl.Phys. B165, 511, 

[31] R. Jackiw and A. Strominger, Phys.Lett. B99, 133 (1981), 

[32] A. di Sessa, J.Math.Phys. 15, 1982 (1974), 

[33] B. Svaiter and N. Svaiter, Revista Brasileira de Fisca 20, 15 (1990). 

[34] A. Makhlin, Phys.Rev. C63, 044903 (2001), arXiv:hep-ph/0007301 [hep-ph]. 


[35] K. Dusling, T. Epelbaum, F. Gelis, and R. Venugopalan, 

arXiv: 1009.4363 [hep-ph]. 

Nucl.Phys. A850, 69 (2011) 

[36] M. Luzum and P. Romatschke, Phys.Rev. C78, 034915 (2008), 

arXiv:0804.4015 [nucl-th]. 


[37] P. Romatschke, Int.J.Mod.Phys. E19, 1 (2010), arXiv:0902.3663 [hep-ph]. 

[38] A. Krasnitz and R. Venugopalan, 

arXiv:hep-ph/9809433 [hep-ph]. 

Nucl.Phys. B557, 237 (1999) 

[39] T. Epelbaum and E. Gelis, Phys.Rev.Lett. Ill, 232301 (2013), 

arXiv: 1307.2214 [hep-ph]. 


[40] 1. Balitsky and G. A. Chirilli, Phys.Rev. D88, 111501 (2013), c 

irXiv: 1309.7644 [hep-ph]. 

[41] A. Kovner, M. Lublinsky, and Y. Mulian, Phys.Rev. D89, 061704 (2014) 


arXiv: 1310.0378 [hep-ph], 

[42] A. Kovner, M. Lublinsky, and Y. Mulian, JHEP 1408, 114 (2014), arXiv:1405.0418 [hep-ph], 

[43] S. Caron-Huot, JHEP 1505, 093 (2015), arXiv:1309.6521 [hep-th], 

[44] S. Caron-Huot, (2015), arXiv:1501.03754 [hep-ph]. 









